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"q ' served in numerical studies of gravitational collapse. These studies suggest, 
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among other things, that black holes initially form with infinitesimal mass. 
We show generally how a black hole mass formula can be extracted from a 



^ ' transcendental equation. 



Using our approach, we give an explicit one parameter set of metrics 
that are asymptotically flat and describe the collapse of apriori unspecified 
but physical matter fields. The black hole mass formula obtained from this 
metric exhibits a mass gap - that is, at the onset of black hole formation, the 
mass is finite and non-zero. 
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A fundamental problem in general relativity is the investigation of the gravitational 
collapse of matter fields. The main motivation for studying this is the cosmic censorship 
conjecture, one form of which states that gravitational collapse always results in a black 
hole. 

Recently the collapse problem has been studied numerically and the results are intriguing. 
For the spherically symmetric collapse of a scalar field, Choptuik JlJ has found that when 
the initial matter field is an ingoing pulse, the collapsing matter forms a black hole with 
mass given by M = K(c — c*) 7 , where K is a constant, c is any one of the parameters in 
the initial data for the matter field, c* is the critical value of this parameter that gives a 
zero mass black hole, and 7 ~ .37. In particular, no black hole is formed when c < c*. An 
important feature of this result is that it appears to be independent of spherical symmetry 
and the type of matter fields: the same mass formula for the black hole has been found 
for the axisymmetric collapse of gravitational waves by Abraham and Evans |J and for the 
spherically symmetric collapse of radiation by Coleman and Evans |J|. 

Thus, these results appear to reflect a universal property of the Einstein equations in 
strong field regions. So far there is no analytical understanding of this result, nor is there 
an explicit metric that exhibits this behaviour. 

There have however been a number of attempts at explanations, most of which involve 
self-similar solutions f|J|. However, the particular solutions discussed do not appear to 
be relevant for the collapse problem with asymptotically flat boundary conditions, because 
they are cosmological, and the collapse results in infinite mass black holes; the radius of 
the apparent horizon tends to infinity with time. Another (non-self-similar) time dependent 
solution, given by the authors, has the same shortcoming ||. The problem has been discussed 
in two dimensions where an exponent of .5 was reported J/]. There have also been proposals 
for using perturbations of black holes to calculate the exponent analytically for the super- 



critical case c > c*, for which the black hole always forms P-Pl- Recently, further numerical 
perspectives on this problem have been given 



In this Letter we present an alternative analytical approach to this problem. The first 



part of the paper gives a general approach that may be used to extract a black hole mass 
formula from a spherically symmetric time dependent metric. The mass formula arises from 
the solutions of a transcendental equation. The second part of the paper gives an explicit 
example of the procedure: We give a metric that satisfies the dominant energy condition 
and describes a realistic collapse. This specific example suggests that at the onset of black 
hole formation, there may be a mass gap. We comment on how our approach may be used 
to get a black hole mass formula without a mass gap. 

We start with the general spherically symmetric line element 

ds 2 = -fe-^dv 2 + 2e^dvdr + r 2 dQ 2 (1) 

where f(r,v) and tp(r,v) are functions of the radial coordinate r, (0 < r < oo), and the 
advanced time coordinate v, (0 < v < oo), and dfl 2 is the line element of the unit two- 
sphere. The parametrization (JJ) has been used to study the collapse of null shells [|nj . The 
mass function m(r, v ) is defined by 

f{r,v) = {l-2m{r,v)/r), (2) 

and is a measure of the mass contained within radius r. The Schwarzschild solution results 
from setting m(r,v) = constant and ip(r,v) = 0. The mass function will in general be 
parametrized by initial data parameters q. For simplicity we will assume dependence on 
only one parameter, c. 

The apparent horizon is the three-surface that separates regions containing trapped sur- 
faces from the normal regions of the spacetime. This surface is defined by the equation 
g ab d a rdbr = 0. For the metric ([!]) it is given by 

T 

m{r,v;c) = ~. (3) 

We are interested in the asymptotic v — » oo solutions of this equation because this gives 
the large time behaviour of the apparent horizon, and this is what has been investigated 
numerically Taking this limit results in the one parameter family of transcendental 

equations 
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T 

M(r;c) := lim m(r,v;c) = — . (4) 

The solutions of this equation give the radial coordinate of the asymptotic (i> — > oo) apparent 
horizon as a function of c. In general there will be ranges of c for which there are no solutions 
to (U). This will give the subcritical region where no black hole forms. Similarly there will 
be critical and supercritical solutions. 

In the numerical work cited above, the black hole mass is defined by the radius of this 
asymptotic apparent horizon. This is a reasonable definition because it is in this limit that 



the apparent horizon approaches the event horizon [15|. In our discussion, this mass is given 



explicitly by the solutions r^fffc) of these transcendental equations, namely 

M BH (c) := — - — . (5) 

A plot of the solution of the equations (f|) as a function of c will give M BH (c). We emphasize 
that the steps outlined above are general in the sense that if one is given an exact collapse 
solution, this procedure may be used to see if critical behaviour exists in the long time limit. 
Our purpose in this paper is to study what specific mass formulas can be obtained from (|J). 

The mass function in the metric @ must satisfy certain physical conditions in order 
to define a realistic collapse. These are: (i) The metric should be asymptotically flat, (ii) 
the mass function should satisfy m > and dm/dr > 0, and should give, in the v —>■ oo 
limit, a set of transcendental equations from which we can obtain Mbh{c), and (iii) the mass 
function should increase with advanced time corresponding to an implosion of matter, at 
least initially. Also, if the mass function leads to multiple apparent horizons, the outermost 
one will serve to extract a mass formula. 

We are interested mainly in the supercritical case, where a black hole always forms at 
the end point of collapse, since the main goal of our work is to extract a black hole mass 
formula. However, as we will see, the subcritical case (no black hole) also arises naturally. 

The Einstein equations G a b = 87iT a b for the metric ([!]) are 

m' = -4vrr 2 T/, (6) 



m = Arrr 2 T v r , (7) 
ip' = A7crT rr , (8) 

where the prime and dot denote r and v derivatives. We will not fix the stress-energy tensor 
by specifying any specific type of matter, but will instead determine this tensor in terms of 
m and ip. The goal then will be to see if these functions can be fixed so that the energy 
conditions for realistic matter are satisfied. (We note that ip is not fixed by the Einstein 
equations). 

As an example of the above procedure, we now construct one metric that describes a re- 
alistic collapse. This requires a mass function satisfying the above three physical conditions, 
with xjj still arbitrary. A suitable choice is 

m(r, v;c)=c [tanh(e c In r) + 1] tanhv . (9) 

The Arnowit-Deser-Misner (ADM) mass for this choice is 2c. 

The one parameter (c) set of transcendental equations (|j) resulting from this in the 
v — *■ oo limit are 

c [tanh(e c lnr) + 1] = ^, (10) 

which give, using (||), the equation for the black hole mass M BH : 

c (tanh[e c \n(2M BH )} + 1) = M BH . (11) 

Equation (|IID may be solved numerically, and there are solutions for c > c* = .465727. 
In the supercritical region c > c*, the solutions are fit approximately by the equation 

M BH = 0.64 + 20.69(c - c*) - 2205.15(c - c„) 2 + 113620(c - c*) 3 . (12) 

A plot of M BB vs. (c — c*) is given in Fig. 1. The solutions of the transcendental equations 
( PD and the fit ( |T2| ) were found using Mathematica, and also checked with Maple. (The 
Mathematica commands we used for Fig. 1 are: tt := Table[{ (c — .465727), (M)/. Find- 
Root[ c (Tanh[ Exp[c] Log[2M]] +1) == M, {M, 0}] ]}, {c, .465728,-475000, .0001}] Fit[ tt, 

{1, X, X 2 , X 3 }, X ].) 
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There are several comments in order for the fit fll2|) : (i) The mass formula is not a 
pure power law as observed in the numerical integrations of the scalar field and null fluid 
]1].|3|]. (ii) This particular example exhibits a mass gap, because at criticality Eqn. ( |TT ) gives 



Mbh = -64, as Fig. 1 clearly shows, (iii) For c < c* there are no solutions of (0), and hence 
no black hole formation, (iv) For c > c* there are two intersections of the mass function with 
the r/2 curve, and hence an inner and outer apparent horizon. The above fit corresponds 
to the outer horizon, which is what is relevant for determining black hole masses. 

In summary, although this mass function example does not give the mass formula ob- 
tained in the full numerical integrations f^-^j it does give the subcritical, critical and su- 
percritical regions. Furthermore, as we will see below, this mass function corresponds to 
matter satisfying the dominant energy conditions, and thus gives a physical asymptotically 
flat collapse solution. This suggests that there are physical solutions that exhibit a mass 
gap. 

So far the metric function tp(r,v), is arbitrary. The task now is to see how this function 
is restricted so that the metric arises from realistic matter. We would like to give at least 
one example of such a ip(r,v) in order to have an explicit metric. To do this we first turn 
to a discussion of the energy conditions. 

The energy conditions may be imposed most simply by first diagonalizing the stress- 
energy tensor, (determined in terms of m and ip), by solving T a b Vb = Xv a . Diagonalizing the 
(f,r) part of T a b gives the eigenvalues A±: 



4 „. A± = _ m - + *M { t ± ^ + ^1 , . (13) 



The corresponding (v , r) components of the eigenvectors are 



•Ml,^), (M) 



m 



and their norms are 



|« (±) | = ( 47rr2A±W )[ 2e^ + /( 4?rr2A± + m/ ) ]. (15) 



m m 
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(The (9, 0) parts of the stress-energy tensor are determined by Gee and = sin 2 9Gee, 
and are already diagonal (since Ge$ = 0). The corresponding eigenvectors are spacelike.) 

Stress-energy tensors are classified by their eigenvectors, and for physical fields this tensor 
must be either Type I, for which there is one timelike and three spacelike eigenvectors, or 
Type II, for which there are two null and two spacelike eigenvectors JT5] . All physical fields 
are of type I, except for certain null fluid flows, which are of type II. 

We will focus on the type I tensors. Let — p be the eigenvalue corresponding to the 
timelike eigenvector, and 7Tj (i = 1, 2, 3) the eigenvalues corresponding to the three spacelike 
eigenvectors. Then the weak energy condition, which requires that the energy density be 
non-negative, is 

p > 0, and iii > — p for i = 1,2, 3. (16) 

The dominant energy condition, which requires that energy flows are never spacelike, imposes 
in addition to (H^), the condition 



Hi < p for i = 1,2,3. (17) 

From our chosen mass function (|9|), we see that rh > 0, which requires that if)' > 
for the square root in the eigenvalues (|13|) to be real. We note from (|13|) and flT5|) that 
Ame~^ / r f 2 tp' > 0, (which is always true), implies that the eigenvector is timelike as 
required. Also, «W is always spacelike. Therefore, for the weak energy condition we need 
a if} that satisfies ip' > 0. In addition, to preserve the asymptotic flatness of the metric we 
require 

lim ib(r,v = oo) = constant. (18) 

r^oo 

Since we already have A + > A_, for the weak energy condition we require, in addition, a if> 
such that 

A- < 0, (19) 

and 



G ee = (r + m)i>' - 3rmV' + r 2 f{ip') 2 - r 2 e~V - rm" + r 2 f<ip" > A_. 



(20) 



For the dominant energy condition we also require 



- A_ > A 



+ 



and — A_ > G, 



(21) 



that is, all the pressures must be bounded between p and —p. 

A if) that satisfies both the weak and dominant energy conditions is of the form 



where the functions / and g are everywhere positive and have a lower bound > 2. This last 
condition is necessary to enforce the dominant energy condition, which was checked to be 
true numerically for r > horizon radius, and for subcritical and supercritical ranges of c. 

Our main result is a general method for obtaining a mass formula for black holes via 
solutions to a transcendental equation: given any explicit collapse solution, a black hole 
mass formula must arise in this way, (and it may or may not be a pure power law). As 
an explicit example of the procedure, we have given a metric that describes a physically 
realistic spherically symmetric collapse, and that exhibits a mass gap at the onset of black 
hole formation. 

It would be of interest to see if there are other choices for the mass function that describe 
a realistic collapse. In particular, it would be of much interest to find the mass functions 
that give gapless black hole mass formulas. 

An important question 0] concerning critical behaviour in our approach is what features 
of the mass function are responsible for giving criticality. The question actually has at least 
two parts: What features give the three subcritical, critical, and supercritical regimes, and 
what features give rise to a pure power law, and hence a critical exponent? The answer to 
the first question is that any mass function that grows with the parameter in approximately 
the same way as ours is sufficient. We speculate that the answer to the second question 
requires an approximately step shaped mass function (like ours) such that its knee remains 
just touching the r/2 curve (Fig. 2) for a wide range of parameter values c, (instead of 




(22) 
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crossing it, and thereby giving an inner horizon also, as ours does). In our approach the 
former choice of mass function, (that just touches the r/2 curve), requires 'fine tuning' 
whereas the second choice appears to be more generic. A deeper question is what properties 
of the Einstein equations in this context give rise to mass functions having these features. 

A further point to note is that critical behaviour may occur even for unrealistic matter: 
we simply leave ip arbitrary, up to requiring asymptotic flatness of the metric, and not worry 
about imposing the energy conditions. 

It seems possible to use this approach for axially-symmetric metrics as well; one could 
for example use the mass function given here, and still have two remaining metric functions 
at hand for satisfying the dominant energy condition. 

We would like to thank Werner Israel, Ted Jacobson, Karel Kuchar, Don Page, Richard 
Price, and Lee Smolin for very helpful comments. The work of V. H. and E. M. was 
supported by the Natural Science and Engineering Research Council of Canada, and that 
of D. N. partly by DGAPA, National A. University of Mexico. 
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FIGURE CAPTIONS 

Figure 1. M BH vs. (c — c*) curve whose best fit is Eqn. (12). 

Figure 2. The critical mass function (c = .46573) together with a typical subcritical case 
(c=.3). The other curve is the function r/2 which touches the critical mass function. 
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